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ABSTRACT 

The high frequency component in blazars is thought to be due to inverse 
Compton scattered radiation. Recent observations by Fermi-LAT are used to 
evaluate the details of the scattering process. A comparison is made between the 
usually assumed single scattering scenario and one in which multiple scatterings 
are energetically important. In the latter case, most of the radiation is emitted 
in the Klein-Nishina limit. It is argued that several of the observed correlations 
defining the blazar sequence are most easily understood in a multiple scattering 
scenario. Observations indicate also that, in such a scenario, the blazar sequence 
is primarily governed by the energy density of relativistic electrons rather than 
that of the seed photons. The pronounced X-ray minimum in the spectral energy 
distribution often observed in the most luminous blazars is discussed. It is shown 
how this feature can be accounted for in a multiple scattering scenario by an 
extension of standard one-zone models. 

Subject headings: BL Lacertae: general — galaxies: active — gamma rays: the- 
ory — radiation mechanisms: non-thermal — scattering 



Introduction 



The spectral energy distribution in blazars exhibits two distinct components. It was rec- 
ognized early on that the low frequency component is almost certainly due to synchrotron 
radiation. The high frequency component is likely to be inverse Compton (IC) scattered ra- 



diation, although various h adronic origins have also be en considered (e.g.. lMannheimlll993 



Miicke &: Protherodl200l[ ). iPadovani fc Giommi (|l995h suggested a classification of blazars 
based on the frequency where the synchrotron component peaks. The distribution of this 
peak frequency showed a rough bimodal distribution, which gave rise to three classes; two 
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classes corresponding to the two largest groups of blazars with peak frequencies in the in- 
frared and in the UV/X-ray, respectively, and a third, smaller class encompassing those 
with intermediate peak frequencies. This rough classification scheme still remains useful 
although the detailed criteria for each class , as well as the nom enclature, have changed. In 
this paper the nomenclature introduced by lAbdo et al.l ( ]2009bl ) is followed. The Low Syn- 
chrotron Peaked blazar (LSP) class includes both flat spectrum radio quasars (FSRQ) and 
low frequency peaked BL Lac objects, while the other two classes, Intermediate Synchrotron 
Peaked blazars (ISP) and High Synchrotron Peaked blazars (HSP), correspond roughly to 
intermediate and high frequency peaked BL Lac objects. 

With the advent of more observational information about the hi gh frequency c o mpon ent, 
this classification scheme was developed into a blazar sequence by iFossati et al.l ( 119981 ). It 
was found that the properties of the IC-component followed a trend reflecting those of the 
synchrotron component; a lower peak frequency of the synchrotron component corresponds 
to a higher luminosity and a lower peak frequency of the IC-component. This could be 
understood as a result of increased cooling in the more luminous blazars. The more extensive 
and detailed observations by Fermi-LAT have shown that these trends not only connect the 
three classes of blazars but actually constitute smooth and continuous relations between 
the properties of th e IC-component and the peak frequency of the synchrotron component 
( lAbdo et al.ll2010bl ). Furthermore, Fermi-LAT has revealed an additional correlation, which 
shows that more luminous blazars tend to have steeper spectra in the LAT frequency range 
( lAbdo et al.ll2009aL lioiOah . 



It is usually agreed that the electrons producing the IC-componen t are responsible also 



for at least the high frequency part of the synchrotron component (e.g..lMaraschi et al.lll992 



Dermer fc Schlickeiserlll993t ISikora et al.lll994j ; iBloom fc Marscherlll996l ). The origin of the 
seed photons for the IC-component is somewhat more uncertain; they could either be external 
photons or the intrinsically produced synchrotron photons (SSC-models). Since there is a 
rough correlatio n between emission lin e strength and luminosity of the IC-component, it has 
been suggested ( IGhisellini et al.l Il998[ ) that external photons dominate the cooling in the 
most luminous blazars, while the least luminous blazars correspond to the SSC-case. It is 
normally assum ed that the dominant contribution comes from fi rst order inverse Compton 
scattering (e.g., IGhisellini fc Tavecchiol 120091 ; ISikora et al.l 120091 ) . The spectral properties 



are then determined mainly by the energy distribution of the injected relativistic electrons. 
Another possibility is that the source properties are such that the IC-component is due to 
multiple inverse Compton scatterings (MIC). In this case, its spectral properties depend 
to a large extent also on the scattering process itself. In addition to different spectral 
characteristics, time variations are expected to have properties which differ between MIC- 
models and first order scattering models; this is particular true for correlations between flux 
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variations in the synchrotron and IC components. 



In iBjornsson fc Aslaksenl (120001 . henceforth BA00) MIC-models were discussed in con- 
nection with the observed properties of the radio/mm outbursts in blazars. Motivated by 
the Fermi-LAT observations of blazars, the present paper extends this discussion to the high 
frequency radiation. Guided by the numerical results in BA00, an analytical description of 
MIC-models is developed in §2. It is emphasized that the origin of the seed photons only 
marginally affects the results. The consequences for the synchrotron component is briefly dis- 
cussed in § 3, in particular the effects of cooling on the synchrotron self-absorption frequency. 
The most luminous blazars often exhibit a distinct X-ray minimum in t heir spectral energy 



distri bution, which is hard to account for in one-zone MIC-models (e.g., iBloom fc Marscher 



19961 ). The extensions of standard one- zone models needed to accomodate such minima in 
MIC-models are discussed in §4. In MIC-models, the relativistic electrons radiate most of 
their energy in the Klein-Nishina limit. Hence, with a typical bulk Lorentz factor T ~ 10 for 
the emission region in blazars, the photons detected by Fermi-LAT shou l d be u p-scattered 
from a freque ncy range overlapp ing that of Swift-BAT. ISambruna et al.l (120101 . henceforth 
S09, see also lAbdo et al.l l2009bl ) have discussed the properties of a sample of blazars de- 
tected both by Fermi-LAT and Swift-BAT. This is taken as the starting point in § 5 for 
an evaluation of how MIC-models and first order scattering models fare in a comparison 
with observations. It is argued that several observed properties are more straightforwardly 
accounted for by a multiple scattering scenario than by one where only a single scattering 
occurs. This has implications for the origin of the blazar sequence, since in a MIC-scenario 
observations indicate that this sequence is not primarily due to variations in the energy 
density of seed photons but rather to the energy density of relativistic electrons. The main 
conclusions are collected in 5 6. 



2. The inverse Compton spectrum 



Zel'dovich (see lRybicki fc Lightmanlll979l ; iPozdnyakov et al.lll976l ) originally proposed 



that the spectrum resulting from multiple inverse Compton scatterings by thermal/mono- 
energetic electrons in the Thomson limit can be approximated globally by a power-law. 
The usefulness of this approximation for describing the multiply inverse Compton scattered 
radiation depends mainly on the number of orders that can be accomodated before the 
Klein-Nishina limit sets in. This in turn is set by a combination of the frequency of the seed 
photons and the energy of the electrons, which give the dominant contribution to the inverse 
Compton flux. For a power-law distribution of relativistic electrons the latter energy needs 
to be determined in a self-consistent way. 
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2.1. Contributions to the inverse Compton flux 

Consider a power-law distribution of flux (f(v)) with spectral index a (f(v) oc v~ a ) 
due to multiple inverse Compton scatterings by an electron distribution ^(7) = fCy _p with 
p > 2 and 7 > 7 m i n , where 7 is the Lorentz factor of the relativistic electrons. Assume further 
that a < 1 so that cooling occurs in the Klein-Nishina limit, i.e., electrons cool on photons 
with frequency v oc 7 . Let the electrons be uniformly injected into a region of size r Q from 
which they escape with speed v esc . Define j c such that its cooling time t c (7 = 7 C ) = r /v csc . 
The cooling time for electrons with 7 > 7 C is then given by 



*c(7 > 7c) / 7c 



*c(7c) V 7 



(1) 



In a steady state situation, the main contribution of photons to the energy density 
U p h{v){oc uf(u)) at frequency v comes from electrons with 7 determined by (cf. BAOO) 



U p h{^i) ( 7\ 3 V Al) ( 7 



£V(zA7c) \7c/ r Q \7 



( 7c 

7 



p-l-2a 

(2) 



The factor oc 7 3 ~ p corresponds to the energy gain (7 2 ) per scattering times the density of 
electrons with Lorentz factor 7 (oc 7 1_p ), while the factor oc 7~ 2 ( 1-a ) takes into account the 
fact the electrons with different Lorentz factors scatter on different photons [y oc 7 -2 ) and, 
hence, different energy densities. The cooling distance r{pj) < r Q for 7 > j c . For 7 C > 7 m i n , 
there are two cases. 

a > (p — l)/2: If 7 > 7 C then r(j)/r = ( r y c /l) a and the ratio in equation (J2]) is 
(7c/7) p_1 " a which is smaller than unity, since p > 2. Hence, electrons with 7 > 7 C cannot be 
the main contributors to the flux in the inverse Compton component. For 7 < 7 C , r(7) = r Q 
and the ratio in equation (T5]) is given by (7 c /7) p-1 ~ 2a , which is also smaller than unity, since 
a > (p — l)/2. Therefore, the flux in the inverse Compton component is dominated by 
electrons with 7 = j c . 

a < (p — l)/2: For 7 > 7 C the result is the same as above. However, when a < (p — l)/2 
the ratio in equation fl2]) is larger than unity for 7 < 7 C . Therefore, in this case, the flux 
in the inverse Compton component is dominated by the electrons with lowest energies, i.e., 
7 = 7min- This case also includes p > 3, since a < 1 is a requirement for the present 
discussion. 

When 7 C < 7 m i n , cooling is important for all the electrons. It is straightforward to 
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show with arguments similar to those above that in this case the inverse Compton flux is 
dominated by electrons with 7 = 7 min independent of the value of a. 



2.2. Cooling characteristics 

The results in § 12.11 indicate that the spectral shape is intimately connected to the 
cooling characteristics of the source, i.e., the value of 7 C . The cooling time for an electron 
with Lorentz factor 7 is 

3mc 

where or is the Thomson cross-section. Since cooling occurs in the Klein-Nishina limit, the 
frequency to use in equation ([3]) is v oc 7 _1 . The cooling characteristics differ somewhat in 
the three cases discussed above. 

a > (p — l)/2: t/ph(^) is determined by electrons with 7 = j c . It can be expressed as 

u P M = [ — ) £W ( 4 ) 



/A, 



where z/ max = / y c mc 2 /h and {7 p h, is the flux resulting from the complete cooling of electrons 
with Lorentz factor 7 C , 

— = f esc e e (7c)f4- (5) 

Here, U e is the total energy density of electrons injected into the source and e e (7 c ) = 
(7c/7min) 2 ~ p is the fraction thereof that has energies around 7 = j c . 

Since t c (7c) = r o/ v esc, equation (j3J) implies 

To = ll~ l - 2 \ (6) 

which relates the values of a and 7 C . Here, 

-7 mi n^e^Tr , (7) 



3(p - 2) 



where n e is the total number density of electrons. Hence, r Q is roughly the Thomson scatter- 
ing depth (r T ) for an electron distribution with 7 min = 1 ( i.e., r D « TTT^Tn ) • F° r consistency, 
this implies r Q < 1, since a > (p — l)/2. 

a < (p — l)/2: f/ p h(^) is determined by electrons with 7 = 7 m i n . The difference is now 
that i/ max = 7 min mc 2 //i and that e e (7 c ) is substituted by (7mm/7c) a , which is the fraction of 
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the energy that electrons with 7 = 7 m ; n lose before escaping. In analogy with the analysis 
above, this leads to 

T o = 7^"n" 2a - ( 8 ) 

Again, consistency implies r D > 1, since a < (p — l)/2. Hence, if r D = 1 occurs for 7 C > 
7mm, ol = (p — l)/2 at this instance. Furthermore, as can be seen from the definition of 
r (cf. eq. [7]), equation (jSJ) is equivalent to r T « 7~f", which is the result originally 
derived by Zel'dovich. The reason is that in this limit, electrons with 7 > 7 min contribute 
only marginally to the inverse Compton spectrum and, hence, the electron distribution 
corresponds effectively to a monoenergetic distribution. This also accounts for the fact that 
a is independent of j c . 

When 7 C < 7 m ; n , U p h(is) is determined by electrons with 7 = 7 m ; n . However, in contrast 
with the previous case, electrons with 7 = 7 min lose all the energy and, hence, the factor 
(7min/7c) a is substituted by unity, which leads to 



Physically, the factor (7 m in/7c) a accounts for the decrease in the column density of electrons 
with 7 = 7 min caused by cooling. 

The value of a describes the upper envelope connecting the different Compton orders 
and is determined by electrons of a given energy (i.e., 7 min or 7 C ) . Hence, the local spectral 
index can deviate substantially from a, in particular, for large values of 7 m ; n or j c . However, 
as discussed in BAOO, a power-law distribution of electron energies considerably smoothens 
the spectrum as compared to a thermal/mono-energetic distribution. A distinct feature of a 
power-law distribution is that the energy of the electrons determining the value of a changes 
rapidly from 7 C to 7 m i n (or vice versa) as the value of r Q goes through unity. Although the 
value of a changes smoothly in this transition, the local spectral index can change rapidly, 
in particular for cases with 7 C ^> 7 min at this instance. This change in spectral character 
from a "bumpy" to a smooth spectrum as 7 C is replaced by 7 min as the main contributor to 
the inverse Compton emission is evident in the spectra presented in BAOO. 

These expressions for a and 7 C for a given value of r Q depend only on the assumption 
of cooling in the Klein- Nishina limit. In order to relate the values for a and 7 C to the 
source properties, the origin of the seed photons needs to be considered. Two different 
scenarios are usually envisaged; external photons, whose properties are unrelated to those 
of the relativistic electrons, and synchrotron photons produced by the relativistic electrons 
themselves (SSC- models). 
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2.3. The seed photons 



When the seed photons are due to synchrotron radiation from the same electrons pro- 
ducing the inverse Compton flux, a particular simple description is possible. This is due to 
the fact that synchrotron radiation can be regarded as inverse Compton scattering of the 
virtual photons associated with the magnetic field. Hence, the second equation needed to 
obtain separate values for a and 7 C is straightforward to derive. 

a > (p— 1)/2: The expression relating the synchrotron and inverse Compton components 

is 

where U p h, is given by equation (j5J). The energy density of the magnetic field is denoted by 
Ub and the cyclotron frequency by v B . This yields 



o \ I— a / \ 2— p 

7 c mc z \ / 7c \ e e 



hVn J \7min/ CB 

where e e /e B = 4:V esc U c /cU B . 

In order to single out the density dependence, it is useful to write 

— = e e BoTo, (12) 
£b 

so that e e Bo is the value of e e /e B for r Q = 1. Convenient expressions for a and 7 C are obtained 
by evaluating equation ffTTl) for r D = 1 (i.e., a = (p — l)/2), which yields 

x (p-l)/2 / 2X (p-3)/2 

/co \ / /mm'"' 1 ^ 



j i t i CeBo, (13) 

,7min/ V hv B J 

This alternative parameter (7 co ) is the value of 7 C corresponding to r Q = 1. Hence, the 
synchrotron-inverse Compton spectrum is basically described by two parameters: The value 
of t , which depends on the column density of electrons but not the magnetic field and 7 co , 
whose value includes the magnetic field but not the column density of electrons. In the 
external photon scenario, the same discussion applies, except that the frequency and energy 
density of the seed photons are not a priori related and, hence, the value of 7 co depends on 
two parameters instead of just one as in the SSC-case. 

With the use of equations ([HD and (TTU|) . separate expressions for a and 7 C can be derived 
in the form of two second order equations (cf. BAOO). However, in order to highlight a few 
characteristics of j c , it is more convenient to express them as 
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and 




In y 
In 7c 



(15) 



where y = mc 2 jhv^ has been introduced. Since for blazars the value of y is likely to be 
very large (e.g., y ~ 10 14 for B ~ 1), it is seen from equation (fT5l) that the value for 7 C is 
expected to be quite sensitive to that of r D . Furthermore, as the value of r D increases, the 
decrease of 7 C accelerates, due to the 1/ hi7 c factor in the exponent. Both of these properties 
derive from the fact that the effective number of inverse Compton orders is proportional to 
In y/ In 7 C . 

In the external photon scenario, equations (TH1) and ffTBI) completely describe the relation 
between the source properties and the values of a and j c . This is not quite the case for the 
SSC-scenario, since it was tacitly assumed in equation (TlO|) that 7 C > 7 a b s , where 7 a b s is 
the Lorentz factor corresponding to the synchrotron self-absorption frequency (f a bs)- When 
7c < 7abs, the energy density of synchrotron photons peaks at i/ abs rather than the cooling 
frequency {y c = 7cV B ). Hence, the same analysis as above can be made except that two 
substitutions need to be made: U-q —> ^ P h( I/ abs) and v-q ^abs = 7abs z/ B- Since 



where the last factor on the RHS accounts for the decreasing column density of electrons 
with 7 = 7abs due to cooling, the expression corresponding to equation (fT5l) is 



When r = 1 occurs for 7 C < 7 a b s , 7abso is the value of 7 a b s at r Q = 1. In this case, separate 
values for 7 co and 7 a bso cannot be obtained; however, their combination, which occurs in 
equation ( 1T7|) . is given by an expression similar to that in equation ( TT5j) 



The main difference between equation (fl5l) and equation (TTTI) is that the exponent of the 
RHS-side is roughly a factor of two smaller for the latter as compared to the former. Formally 
the reason for this is that when r Q = 1 occurs for 7 C > 7 abs , 7 a bso — ^ 7co and equation ffT8]) 
becomes identical to equation ffl3|) . Hence, the variation of 7 C with r Q slows down considerably 
for 7 C < 7ab s . The physical reason for this is that the energy density of synchrotron seed 
photons varies much slower in this regime. This conclusion hinges on the assumption that 




(16) 




(17) 




(18) 
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the value of 7 a b s varies more slowly than j c . It is shown in the next section that this is 
indeed the case. 

a < (p — 1)/2: In this regime, a is obtained directly from equation (jSJ). The variation of 
7 C is quite similar to that for a > (p — l)/2 and expressions corresponding to equations (IT51) 
and (ITTj) can be derived. The main difference is that 7 C is replaced by 7 m i n in the exponent 
of t , which is due to the fact that the main contribution to the inverse Compton flux now 
comes from electrons with 7 = 7 m i n . These results also remain valid for the regime 7 C < 7 m i n - 



3. The synchrotron spectrum 

In order to relate the properties of the synchrotron component to those of the inverse 
Compton component, the value of 7 a b s needs to be determined self-consistently. Since the 
main parameter in the multiple inverse Compton scenario is r D , in this section attention is 
restricted to the case when only the column density of electrons varies. 



When 7 C > 7 a b s , standard synchrotron theory gives (e.g.. |Pacholczyklll970l ) 



7absOC T™. (19) 

Cooling affects the absorption properties for 7 C < 7 a bs, since the column density of electrons 
with Lorentz factor 7 abs varies as (7 c /7abs) a - Hence, the expression for 7 abs in this case is 
obtained by letting r G — > r Q (■y c /labs) a in equation (fl9l) . which leads to 

7absOC (r o7c Q )rai. (20) 

The expression for 7 C to use in equation (120]) is that given either in equation (Tl5!) for the 
external photon case or equation (TTT|) for the SSC-case. For external seed photons, this 
yields 

(l-lny/2ln-, c ) 

7abs0CT o 4+a+P , (21) 

where equation ([6]) has been used. As mentioned above, for a < (p — l)/2 the expression for 
7abs is quite similar to the one for a > (p — l)/2, except that 7 C is replaced by 7 m i n . In the 
SSC-case, the variation of 7 C with r Q is considerably slower than for external photons (cf. 
eqns. [T3] and [17]) and, hence, the exponent in equation ( 12 ip is correspondingly smaller. 
However, the actual value of the exponent is likely to be rather similar in the two cases, since 
the value of y is expected to be much larger in the SSC compared to the external 

photon case (e.g., y ~ 10 14 for B ~ 1 as compared to y ~ 10 5 for optical seed photons). 
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With 7 C < 7 a bs ~ 10 2 , it is likely that \ny > 2 hi7 c independent of the origin of the seed 
photons. The dependence of 7 a b s on r Q would then reverse as compared to the case 7 C > 7 a b s ; 
instead of having a value of 7 a b s increasing with r D , it now decreases. 



4. Connection to blazars 

In a MIC-scenario, most of the main characteristics of the inverse Compton component 
are quite insensitive to the origin of the seed photons and are instead the result of the 
scattering process itself. An example is the expected spectral correlations between the part 
determined by a (i.e., mainly X-ray) and the spectral range where the emission occurs in 
the Klein- Nishina limit (i.e., mainly gamma-ray). 

These spectral correlations can be expressed in terms of the spectral breaks that occur 
at frequencies above the X-ray regime (see BA00). There will be one or two breaks depending 
on the value of a. For a > (p — l)/2, the X-ray spectrum is determined by 7 C and extends 
to hv ft '~ y c mc 2 . Above this frequenc y, cooling is important and the spectral index is given 



by p — 1 (IBlumenthal &: Gouldlll970l ); hence, the magnitude of this spectral break is Aa = 
p—1 — a. This is also the frequency where the value of vf(v) peaks. For a < (p — l)/2, this 
single break splits into two. The X-ray spectrum is now determined by 7 m i n and extends to 
hv ft 7 min mc 2 . Above this frequency, emission occurs in the Klein-Nishina limit implying 
that the emitted frequency is proportional to electron energy. Since the electrons lose only 
a fraction of their energy (oc 7°), the spectral index is given by p — 1 — a, which corresponds 
to a spectral break of magnitude Aa± = p — 1 — 2a. A second, cooling break at hv ft 7 c mc 2 
has magnitude Aa 2 = a. The peak in vf(v) occurs at the first break when a < p — 2 
and at the second break when a > p — 2; the latter requires p < 3. The two breaks merge 
again into one for 7 C < 7 m i n . These spectral properties are summarized in Figured], which 
shows schematically how the inverse Compton scattered part of the spectrum changes as the 
column density of electrons (i.e., r Q ) increases. It may be noticed that of the three possible 
frequency ranges with different spectral indices only the one at the highest frequencies is 
determined exclusively by the energy distribution of the injected electrons (i.e., p). The 
other two spectral indices can vary continuously without any change in p. Such variations 
are instead caused by changes in the electron density and are correlated through the value 
of a. 

The source properties of blazars deduced from observations in a MIC-scenario differ 
in some crucial aspects from those in a single scattering scenario; however, several remain 
the same. In a MIC-scenario, the value of 7 m i n should be small to assure a smooth inverse 
Compton spectrum; from BA00 this implies 7 min < 10. Similar values are also suggested for 
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single scattering of photons from the broad line region, since the X-ray spectra measured by 
Swift-BAT is o n average somewhat too soft to be due to the cooling tail of higher energy 



electrons (e.g., ISambruna et all l2010l ; lAbdo et al.l l2009b[ ). For the most likely site of the 



emission region in a single scattering scenario, cooling on external photons leads to 7, 



< 



10 2 i n the most luminous blazars (LSPs) (e.g., iGhisellini fc Tavecchiol 120091 ; ISikora et al. 
20091 ). This is also similar to the MIC-case and implies that the peak in the spectral energy 
distribution of the synchrotron component occurs close to the self-absorption frequency. 

The main difference between the two scenarios concerns the value of e e /eB- Single scat- 
tering models are usually consistent with rough equipartition between the energy densities in 
electrons and magnetic field (i.e., e e ~ €b), while multiple scattering models require e e ^> cb 
(cf. eq. [13]). Furthermore, the smooth power-law X-ray spectra in LSPs indicate r Q > 1 
in MIC-models in order for this radiation to be produced by the lowest energy electrons (cf. 
eq. 0). 

The distinct minimum in the spectral energy distribution occurring at the transition be- 
tween the synchrotron and inverse Compton components, which is particularly pronounced 
in LSPs, is not consistent with a strictly homogeneous one-zone MIC-model. The low value 
of 7min needed to produce the observed smooth power-law in the X-ray regime would also 
fill in this minimum by inverse Compton scattering of the synchrotron component (e.g., 
Bloom &: Marscherlll996l ). However, any one-zone synchrotron/ Compton model have prob- 
lems accounting for some of the main char acteristics of LSP; for example, although c orrelated 
time variations are sometimes seen (e.g., iMarscher et al.l 12008c lAbdo et al.l 12010a ) between 
the optical and the gamma frequency ranges, this is not generally the case. 

The similar value of the magnetic field B in single and multiple scattering scenarios have 
two implications: (i) The deduced location of the emission site should be roughly the same, 
(ii) In LSPs most of the radiation is emitted by cooling electrons. Hence, the much larger 
value of the energy density of electrons in a MIC-model as compared to a single scattering 
model necessitates a correspondingly smaller emission volume, i.e., r Q <C R, where R is the 
size of the available emission region. This suggests a possible extension of the one-zone MIC 
model that can account for the minimum in the spectral energy distribution. 

For simplicity, assume that the values of B as well cLS ^esc the same throughout the 
emission region, including the smaller injection region. The uniform injection of relativistic 
electrons inside r Q implies that the fraction of electrons with 7 > j c , which escapes and fills 
the whole emission region R, is proportional to their cooling time. The contributions to the 
optically thin synchrotron flux is determined by the average time spent by an electron in 
the two regions. With no additional cooling outside the injection region, the shape of the 
synchrotron spectrum emitted by the escaping electrons is, therefore, similar to that inside r Q , 
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while the flux is a factor R/r Q larger. Although cooling on the inverse Compton component 
is less efficient by a factor R/r Q outside as compared to inside the injection region, cooling 
may steepen the synchrotron spectrum from the escaping electrons at higher frequencies. 
A further contribution could come from synchrotron cooling. Hence, the smaller injection 
region contributes most of the inverse Compton flux, while the synchrotron emission comes 
mainly from the emission region as a whole. 

Within the injection region, the amplitude between consecutive Compton orders is 
roughly 7mi n ' r T ~ 7min ( c ^- e 1- 0)- Outside the injection region this amplitude is 
smaller by a factor R/r Q due to the smaller column density of electrons. The small value 
of 7 C expected in LSPs causes the inverse Compton scattering also outside the injection 
region to be dominated by the low energy end of the electron distribution. The minimum 
in the spectral energy distribution observed in LSPs is typically a factor 10 below the peak 
of the synchrotron component. In order to make MIC-models compatible with such min- 
ima, 7^ a \ Q /R < 1CT 1 /™ is required, where n is the Compton order corresponding to the 
frequency of the minimum. Since 7 m i n < 10 is needed for a smooth X-ray spectrum, the 
strongest constraint on r Q /R is obtained by using 7 min 10, which implies n ~ 2 (taking 
the minimum to occ ur at a frequency 10 17 Hz and a synchrotro n peak frequency at 10 13 Hz). 



With a « 0.6 (e.g., ISambruna et all 1201(1 lAbdo et alJl2009bh . this shows that R/r Q > 20 



gives rise to an X-ray minimum deep enough to be consistent with observations. 

The requirement for a much larger value of e e /e^ inside the injection region than out- 
side constrains the mechanism whereby the escaping electrons are dispersed throughout the 
synchrotron emission volume. Adiabatic expansion is unlikely since the value of ee as well 
as the radiative efficiency would decrease. A constant value of B and no adiabatic losses of 
the electrons, as in the above example, are consistent with scattering possibly in connection 
with streaming. In addition, the acceleration of electrons could occur in regions with low 
values of B, suggesting that the injected e lectrons may derive th eir energy from magnetic 



reconnection (for example, as discussion by lNalewajko et al.ll2010l ). Both of these processes 



may take place simultaneously and, hence, the transition from the synchrotron to the inverse 
Compton component in the spectral energy distribution is sensitive to the combined effects 
of the acceleration and the subsequent spreading of the electrons through the whole emission 
volume. 

In a MIC-scenario, spectral properties and flux variations are quite different from those 
in a source where the high frequency emission is mainly due to a single scattering. As 
discussed above, the actual source properties need not differ too much between the two cases 
due to the sensitivity of 7 C to the density of relativistic electrons, in particular, but also to the 
magnetic field strength/density of external photons. Furthermore, this sensitivity becomes 
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more pronounced for lower values of 7 C (cf. eq. [15]). Even for a continuous distribution 
of intrinsic source parameters one may then expect a bi-model distribution of observed 
properties, although transition cases should be observed. The first group could be adequately 
described by a standard synchrotron/Compton model, in which multiple inverse Compton 
scatterings may occur but, most importantly, emission in the Klein- Nishina limit would 
not dominate the radiation energetically. The second group corresponds to the situation 
discussed above and is characterized by 7 C < 7 a b s - This latter group would have three 
distinct observational characteristics: (1) The emitted energy in the IC-component should 
dominate that in the synchrotron component. (2) With 7 a b s ~ 10 2 , the peak of vf{y) in the 
IC-component should occur at an energy < 10 2 MeV as measured in the rest frame of the 
source. (3) The optically thin synchrotron flux is emitted by electrons with 7 > j c . Hence, 
the peak of vf{y) in the synchrotron component occurs at v « v^- 

The properties of these two groups are reminiscent of those observed for the HSPs (first 
group) and LSPs (second group). In the simplest version, the main difference between their 
intrinsic source properties could be a higher value of e e /eB in the injection region of LSPs as 
compared to HSPs. In this scenario, flux variations in the two classes of blazars are expected 
to have different properties. Since inverse Compton cooling plays a minor role in HSPs, 
variations in the synchrotron component should correlate with those in the IC-component, 
while in LSPs such correlations are expected to be more complex (see §5). 



5. Fermi blazars 



The LAT experiment (lAtwood et al.l 120091 ) aboard the Fermi Gamma-ray Observatory 



has provided a wealth of data related to the high frequency properties of blazars. With 
a typical bulk Lorentz factor T ~ 10 for the emission region in blazars and emission in 
the Klein-Nishina limit, the photons detected by LAT are those up-scattered from energies 
corresponding roughly to the sensitivity range of the Swift-B AT instrument. Hence, within 
the MIC-scenario, conclusions drawn from a comparison of the observed properties of blazars 
detected both by LAT and BAT should be rather model independent. S09 have discussed the 
spectral and luminosity correlations for such a sample of blazars. Since the LAT and BAT 
observations are not simultaneous, a comparison for individual objects is not warranted. 
However, S09 find that the different subgroups (i.e., LSP, ISP and HSP) have distinctly 
different spectral properties. 

For the LSPs they find an average LAT spectral index «g ~ 1.6 and an average BAT 
spectral index ax ~ 0.6. In the MIC-scenario, LSPs correspond to blazars with 7 C < 7 a b s - It 
is therefore likely that the LAT detected emission is due to cooling electrons, which leads to 
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p = 1 + etc ~ 2.6. Since this implies ax < (p — l)/2, the X-ray spectrum should be due to 
electrons with the lowest energies (i.e., 7 m i n ). Two breaks are expected for 7 min < j c , which 
would result in a rather flat peak in vf(y). The gamma-ray emission from HSPs is usually 
interpreted as first order Compton scattering. Although the number of HSPs in the S09 
sample is low, their mean value of oq (~ 0.8) is consistent with the same value of p as for 
the LSPs and negligible cooling. Furthermore, the spectral properties of the ISPs, which fall 
in between those of LSPs and HSPs, are consistent with being transition cases for which the 
gamma-ray emission is produced in the Klein-Nishina limit but with the cooling frequency 
in or somewhat above the LAT energy range. 

It is well known that in general the p eak frequency of the IC-component shifts to lower 



frequencies for higher luminosities (e.g., iFossati et all Il998l ). S09 show that on a more 



detailed level this is due to more luminous sources having spectra with smaller values of ax 
as well as larger values of chg- Furthermore, this correlation seems to be valid also within 
the group of the most luminous sources, which have their peak frequencies in between the 
BAT and LAT frequency ranges. As can be seen from Figured! this is consistent with a 
MIC-scenario. 

Hence, the varying spectral properties along the blazar sequence can be interpreted as a 
continuous variation o f the cooling f r equen cy. Although this agrees with the suggestion orig- 
inally put forward by IFossati et all (119981 ) . the implications of this in the MIC-scenario are 
quite different. The blazar sequence is normally thought to be governed by the energy density 
of external photons so that a higher density would correspond to a lower cooling frequency 
and, hence, a higher total luminosity. At some luminosity, the energy density of external 
photons have decreased to the extent that int ernally produced synchrotron photons would 
take over as seed photons (e.g., iBottcherl 120071 ) . The former group corresponds to LSPs and 
the latter to HSPs. LAT observations show that LSPs tend to have «g > 1 and ISPs /HSP s 
a G 1) with most of the objects with a G w 1 being classified as ISPs (lAbdo et al.ll2010al ). 
This indicates that the transition from external photons to internal synchrotron photons as 
the dominant source of seed photons corresponds to a cooling frequency lying in or close 
to the LAT frequency range. Such a coincidence remains to be explained in the standard 
interpretation of the blazar sequence but is expected in the MIC-scenario. 

With a bulk Lorentz factor T ~ 10 and emission in the Klein-Nishina limit, the LAT 
spectral range corresponds, roughly, to electrons with Lorentz factors 7 ~ 10 3 . Standard syn- 
chrotron theory implies that synchrotron self- absorption sets in for electrons with 7 a b s ~ 10 2 . 
In the MIC-scenario, the peak of the synchrotron component in LSPs corresponds, roughly, 
to the self-absorption frequency. Hence, electrons radiating in the LAT-frequency range 
have energies similar to those producing the optical synchrotron radiation. The ISPs were 
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originally identified as those blazars having their synchrotron peak frequency in between the 
LSPs and HSPs (i.e., in the optical/near- UV range). Since this peak frequency corresponds 
to 7 ~ 7c, the cooling frequency in the IC-component for the ISPs is expected to fall in the 
LAT energy range. 

Increasing the energy density of the seed photons (i.e., the external photons and/or the 
magnetic field) would decrease the cooling frequency. As can be seen from equations ([6]) and 
dHJ), this would cause the value of a to either increase (a > [p — 1]/2) or remain constant (a < 
[p — l]/2). However, S09 show that «x (~ at) actually decreases with increasing luminosity. 
This implies, in the MIC-scenario, that the main parameter determining the blazar sequence 
is the density of relativistic electrons (i.e., To). Therefore, the salient features of the blazar 
sequence can be reproduced by varying only To (cf. Fig.[T]). The rough correlation observed 
between total luminosity and strength of the emission lines in blazars would then not be the 
primary cause for the blazar sequence but rather a secondary effect induced by a correlation, 
for example, between the jet power and the accretion power (i.e., disk and line emission). 

The spectral indices used by S09 were ca lculated assuming a power-law distribution 



of the flux o yer the whole LAT energy range (lAbdo et al.ll2009al ). However, as shown by 



Abdo et al.l (12010a ). a spectral steepening at the highest frequencies is common for LSPs 
and such a behavior is also seen for some ISPs. Although the steepening can be quite strong, 
it normally occurs at a frequency such that the spectral indices for the lower frequencies are 
only marginally smaller (Aoq < 0.1) than those deduced assuming one power-law for the 
whole spectral range. These small changes in the value of do not significantly change the 
conclusions of the discussion above; however, the origin of the break for the LSPs remains to 
be explained. It may be a cooling break akin to the ones in ISPs, in which case a somewhat 
larger bulk Lorentz factor is implied for LSPs as compared to ISPs. Furthermore, since p is 
now determined by the spectral index pertaining to the frequency range above the break, its 
value will be significantly larger than deduced above and, typically, p > 3 is needed. Such a 
value is too large to account for the non-cooling spectral index for the HSPs and, hence, a 
varying value of p along the blazar sequence is required. 

Another possible cause for the spectral break is absorption on extern al emission line pho- 



tons. This process is normally thought to be dominated by Ly a photons (IGhisellini fc Tavecchio 



20091 ). in which case the expected break occurs at frequencies too high to be compatible with 
observations. However, the corresponding resonance line in He II can be as strong as Ly a in 
photoionization models with spectra as hard as those in quasars. The absorption then sets 
in at a frequency a factor of 4 smaller and would be roughly consistent with the observed 
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breaks in LSPfl The exact break frequency as well as the spectral shape above the break 
would depend on both the geometry of the emission line region and its location relative the 
blazar emission region. 

The correlations discussed so far have assumed steady state conditions. Such an as- 
sumption is not justified during times of rapid changes. Along t he blazar sequence, th e 



value of «g increases with blazar luminosity. It has been noted by iGhisellini et al.l ( 120091 ). 
that such a correlation is not generally valid for variations in individual blazars; in fact, 
observations show that in outbursts there is a tendency for variations in the high frequency 
range to be roughly orthogonal to the one defined by the blazar sequence. Sometimes the 
high frequency spectral peak even moves from below to above the LAT spectral range, i.e., 
the spectral properties of LSPs become more ISPs/HSPs like. If the blazar sequence reflects 
steady state conditions, outbursts could then correspond to variations in the source short 
enough for cooling effects to become less important in the LAT spectral range. 

The much smaller emission region of the inverse Compton component as compared to 
the synchrotron component would suggest a correspondingly smaller variability times scale as 
well as flux variations in the former to precede those in the latter. However, this assumes that 
the variability timescale is long enough for steady state conditions to be reached in the inverse 
Compton emission region. As just discussed, there are indications that rapid outbursts in 
LSPs in a MIC-scenario occur under non-steady state conditions resulting in less severe 
inverse Compton cooling. The fraction of the injected energy emitted as synchrotron emission 
would therefore increase. Furthermore, a substantial part of the synchrotron emission could 
be emitted within the injection region leading to similar variability times scales for the two 
flux components as well as a considerably smaller time lag. Hence, if the correlated time 
variations sometimes seen in the synchrotron and inverse Compton components are due to 
non-steady state conditions, they are expected to be accompanied by changes in the high 
frequency range, which is not along the blazar sequence. Another effect that could smear out 
the expected differences between the variability properties of the two flux components, even 
under steady state conditions, is the structure of the injection region. It has been tacitly 
assumed that the injection of energy occurs in a single region. However, it could take place 
in several regions located throughout the synchrotron emission volume; for example, a quasi- 
homogeneous distribution of injection regions should result in similar variability properties 
for the synchrotron and inverse Compton components. 



1 After the submission of the present paper, a similar suggestion was made by Poutanen & Stern (2010, 
larXiv:1005.3"792ll 
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6. Conclusions 

Multiple scatterings give rise to observed properties quite different from those expected 
in a single scattering scenario. This is mainly due to the different relations between observed 
photon frequency and electron energy; for example, in the case of multiple scatterings, the 
emission over a wide frequency range can be produced by electrons with the same energy, 
while in the single scattering case, there is a direct correlation between emitted frequency 
and electron energy. With respect to the general properties of multiple scattering scenarios, 
the main points of the present paper are: 

• The description of the multiple scattering scenario for a thermal/mono-energetic distribu- 
tion of electrons is extended to also include a power-law distribution. 

• A power-law distribution of electron energies introduces an extra spectral break. 

• The various spectral indices are determined not only by the slope of the power-law dis- 
tribution but also by the density of electrons; hence, spectral differences do not necessarily 
indicate intrinsically different electron distributions. 

• The correlation of time variations in the IC-component with those in a low-frequency 
synchrotron component is expected to change character when cooling becomes important 
for electrons radiating close to or below the synchrotron self-absorption frequency. 

When the observed emission from blazars is interpreted within a multiple scattering 
scenario, the deduced source properties differ in some respects from those derived from a 
single scattering scenario. The main implications from a multiple scattering scenario are: 

• The rate of radiative cooling is more sensitive to the density of seed photons than for single 
scatterings. Even with a rather homogeneous distribution of intrinsic source properties, this 
non-linearity suggests that most blazars would belong to one of two groups, characterized 
either by (i) cooling in the multiple scattering regime is not important or (ii) multiple 
scatterings cause cooling to be important even for electrons radiating close to or below the 
synchrotron self-absorption frequency. This bimodal distribution is reminiscent of the blazar 
classes HSP and LSP, respectively. 

• When multiple scatterings are energetically important, the high energy radiation is emitted 
in the Klein-Nishina limit. Electrons radiating in the Fermi-LAT frequency range would 
then also emit synchrotron radiation in the optical/near-UV range. This implies that blazars 
with their synchrotron cooling frequency in the optical/near-UV should have a corresponding 
cooling break in the Fermi-LAT frequency range. The ISP blazars have properties consistent 
with being such transition objects. 
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• The observed anti-correlation between the total luminosity of the IC-component and the 
spectral index in the X-ray/ Swift-B AT range is best understood as resulting from differences 
in the energy density of electrons rather than that of the seed photons. In fact, it is possible 
that the whole blazar sequence is primarily driven by variations in the energy density of 
electrons. 

• One-zone models are not adequate to describe the spectral transition region between the 
synchrotron and IC components in blazars. Instead a source structure is suggested where in- 
jection of relativistic electrons occurs in a small fraction of the synchrotron emission volume. 
The IC-component is produced mainly within the injection region(s), while the synchrotron 
emission is due to electrons, which have escaped from the injection region(s). Such a source 
structure preserves the main spectral correlations expected between the synchrotron and 
IC-components in one- zone models. 

This research was supported by a grant from the Swedish Natural Science Research 
Council. It was done as part of the International Team collaboration number 160 sponsored 
by the International Space Science Institute (ISSI) in Switzerland. 
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Fig. 1. — Schematic representation of the four qualitatively different energy spectra possible 
in a multiple scattering scenario when a < 1. All the electrons radiate most of their energy 
in the Klein-Nishina limit. The cases 1—4 are distinguished by an increasing value of the 
density of relativistic electrons (r Q ). The spectra are labeled both by their value of 7 C and 
the corresponding range of values for a. Casel: One spectral break; spectral energy peaks 
at the cooling frequency. Case 2: Two spectral breaks; spectral energy peaks at the cooling 
frequency. Case 3: Two spectral breaks; spectral energy peaks at the minimum Klein-Nishina 
frequency. Case 4: One spectral break; spectral energy peaks at the minimum Klein-Nishina 
frequency. For p > 3 only cases 3 and 4 are relevant. 



